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I We have obtained detailed imaging Fabry- Perot observations of the nearby galaxy 

' M82 in order to understand the physical association between the high-velocity 

■ outflow and the starburst nucleus. The high spatial and kinematic resolution of our 

■ observations has allowed us to perform photometric analyses of Ha, [N II], and [O III] 
Q_|. spectral lines at roughly one hundred thousand positions across the extent of the 
Q '. galaxy. 

,H I The observed velocities of the emitting gas in M82 reveal a bipolar outflow of 

^ I material, originating from the bright starburst regions in the galaxy's inner disk, but 

' misaligned with respect to the galaxy spin axis. The deprojected outflow velocity 

indicated by the optical filaments increases with radius from 525 to 655 km s^^. All 
three spectral lines show double components in the centers of the outflowing lobes, 
with the Ha line split by ~300 km s~^ over a region almost a kiloparsec in size. The 
filamentary lobes lie along an axis tilted by 15° with respect to the spin axis, confirmed 
by the regions of line splitting and the ionization pattern over the outfiow. The 
filaments are not simple surfaces of revolution, nor is the emission distributed evenly 
over the surfaces. We model these lobes as a composite of cylindrical and conical 
structures, collimated in the inner ~500 pc but expanding at a larger opening angle 
of ~25° beyond that radius. We compare our kinematic model with simulations of 
starburst-driven winds in which disk material surrounding the source is entrained by 
the wind. There is some evidence for rotation of the wind filaments about the outfiow 
axis in support of entrainment, and we find strong similarities between the observed 
and predicted structures. 
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The data reveal a remarkably low [N II] /Ha ratio in the region of the outflow, 
indicating that photoionization by the nuclear starburst may play a significant role in 
the excitation of the optical filament gas, particularly near the nucleus. An increase in 
the [O III] /Ha ratio along the outflow is observed. At larger radii, the line diagnostics 
and a strong spatial correlation between Ha and soft x-ray filaments are consistent 
with shock ionization. 

A smooth spherical halo around M82 is observed in emission lines, extending to 
at least 2 kpc. We propose that the dusty halo is the primary source of the linearly 
polarized optical emission. A diffuse ionized medium (DIM) with enhanced [N II] /Ha 
emission pervades throughout the stellar disk. We discuss likely sources of ionization 
and heating. 

Subject headings: galaxies: individual (M82, NGC 3034); galaxies: starburst; ISM: jets 
and outflows; galaxies: kinematics and dynamics; techniques: interferometric 



1. Introduction 



Ten years ago Chevalier and Clegg (1985) published a brief letter in which they conjectured 
that a high supernova rate in a galactic nucleus could heat the surrounding gas to temperatures 
with sound speeds exceeding the escape velocity of the galaxy. This hot, tenuous gas would expand 
outward from the galaxy in the form of a "wind," enriching the surrounding intergalactic medium. 
Since that time, the galactic wind model has been confirmed by observations of over a dozen 
galactic-scale outflows: e.g., M82 (|Axon fc Taylor I978|) , NGC 253 ([Fabbiano &: Trinchieri 1981 ), 
NGC 1569 (peckman et al. 1995|) , NGC 1705 QMeurer, Freeman, fc Dopita 1989| ), NGC 1808 
QPhillips 19931) , NGC 2146 (|Armus et al. 1995| ), NGC 3079 ( |Veilleux et al. 19941) , NGC 3628 
([Fabbiano, Heckman, k Keel 1990|) , NGC 4051 ( |Christopoulou et al. 1997| ), NGC 4666 ( pahlem 
et al. 1997] ), Mk 509 (IPhilhps et al. 19831) . These massive ejections of gas and dust are usually 
observed as large (~few kpc), roughly conical structures of filaments originating in the nuclear 
regions and oriented along the minor axes of the galaxies. They are therefore most often observed 
in edge-on disk systems (e.g., M82, NGC 1808), in optical emission lines ( p^ehnert &: Heckman 
1995), radio continuum ( |Baum et al. 19931 ), and soft X -rays (e.g, [Bregman, Schulman, fc Tomisaka] 
1995 ; Armus et al. 1995| ). Recent spectral studies have also shown the ability to detect winds in 



face-on galaxies using their kinematic signatures (e.g., NGC 2782 | Boer, Schulz, &: Keel 1992 ], 



Mk 231 jHamilton fc Keel 1987] ; [Krabbe et al. 19971 ) 



In some cases these winds are indeed driven by supernovae and massive stellar winds from a 
central starburst (e.g., M82, NGC 1569; Lehnert fc Heckman 1995 ; Lehnert fc Heckman 199(^ ), 
while other winds appear to be powered by more exotic forces associated with the central engines 
of AGN (e.g., NGC 3079, NGC 4051; [Colbert et al. 1996a| ; [Colbert et al. 1996b[ ), and some galaxies 
appear to exhibit both starburst and AGN characteristics (e.g., NGC 1808 [[Forbes, Boisson, fc 
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Ward 1992], Mk 231 [ Lipari, Colina, fc Macchetto 1994[ ). Because of the low densities of the wind 



material, emission from optical lines is often very difficult to detect, except in nearby galaxies. 
The x-ray emission from these winds is comparable in luminosity to the optical emission lines, but 
is visible on larger spatial scales and therefore detectable to greater distances. 

Due to its proximity and favorable inclination angle (see Tab. |l]), the irregular disk galaxy 
M82 (NGC 3034) has been studied extensively as a prototype galactic wind system. Following the 
original discovery ( Lynds &: Sandage 1963| ), the first detailed spectroscopic study of the optical 



emission line filaments ( Burbidge, Burbidge, &: Rubin 1964 ) revealed kinematics indicative of a 
bipolar, roughly conical, outfiow of gas along the minor axis of the galaxy. 

Thirty years and over 500 published papers later, the initial interpretation of the M82 
filaments in terms of an outfiow still stands. Support for this picture extends from radio to gamma 
ray wavelengths. The starburst nature of the nucleus of M82 has been verified through its strong 



infrared emission (e.g., Soifer, Houck, &: Neugebauer 1987; Rice et al. 198g; Telesco et al. 1991) 



and numerous compact radio supernovae (e.g., Kronberg, Biermann, & Schwab 1985; Huang et al\ 



1994; Muxlow et al. 1994). Models of the starburst evolution produce appropriate quantities of 



energy and mass on plausible timescales to create and sustain the observed nuclear and galactic 
wind behavior (e.g., Rieke et al. 1993; Poane &: Mathews 1993| ). The wind itself has been observed 
optically in both spectral (e.g., McCarthy, Heckman, fc van Breugel 1987 ; McKeith et al. 1995| ) 
and imaging studies (e.g., Ichikawa et al. 1994| ), in emission lines (e.g., Armus, Heckman, fc Miley| 



1990) and broadband radiation (e.g., O'Connell &: Mangano 1978 ). In particular, kinematic 



evidence for the existence of a galactic wind in M82 has been presented by several optical emission 



line studies (e.g., Heckathorn 1972; Axon & Taylor 197S; Bland & Tully 1988; Heckman, Armus 



fc Miley 1990| ; [McKeith et al. 1995]) , and even by molecular observations ([Nakai et al. 1987| ). 

An extensive x-ray halo has been observed oriented along the minor axis, extending 5-6 kpc 



from the disk, far beyond the visible extent of the optical filaments (e.g., Watson, Stanger, & 



Griffiths 1984; Schaaf et al. 198S; Bregman, Schulman, fc Tomisaka 1995| ). A few percent of the 
supernova energy from the starburst has been deposited in this hot (T ~ 10^ K) "x-ray wind." A 



large (r ~ 8 kpc) spherical halo has been reported at radio continuum wavelengths ( [Seaquist & 
Odegard 1991]) , and has been interpreted as synchrotron emission from relativistic electrons in the 



outfiowing wind. 

With this wealth of data, it is surprising that little attempt has been made to undertake 
detailed comparisons of models and observations. Large-scale galactic winds were first proposed 
to account for the lack of gas in elliptical galaxies ( Mathews fc: Baker 1971 ). The theory of 
these winds has since evolved in tandem with research on starburst-driven galactic winds (e.g.. 
White fc Chevalier 198^ ; ^mith. Kennel, fc Coroniti 1993a| ; ^mith. Kennel, k Coroniti 1993b| ). 
Since the original starburst wind model ( phevalier &: Clegg 1985| ), advances have been made in 
both analytical studies (e.g., Koo fc McKee 1992a| ; [Koo &: McKee 1992b| ) and hydrodynamic 
simulations (e.g., Tomisaka & Ikeuchi 1988; Tomisaka fc Bregman 1993| puchkov et al. 1994 ; 
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Suchkov et al. 1996). Models of winds in other astrophysical situations, such as stellar winds (e.g. 
Shu et al. 1991D and winds from AGN (e.g., [Smith 199"^ ; |Arav fc Li 1994| ; [Arav, Li, fc Begelmaii 



1994), have contributed to our understanding as well. On a broader scale, starburst-driven 



winds have important ramifications for a variety of astrophysical and cosmological situations, 
including enrichment of the intergalactic medium, contributions to the diffuse X-ray background, 
the evolution of dwarf and interacting galaxies, and the formation of elliptical galaxies through 



mergers (see [Heckman, Armus, fc: Miley 1990| and references therein) 



Toward the goal of investigating this model, we have obtained high-resolution imaging 
Fabry- Perot observations of M82 in several optical emission lines. In Section |2| of this paper we 
present our Fabry-Perot observations and describe the methods of data reduction. In Section ^, we 
present the derived two-dimensional emission-line maps of M82 illustrating the distribution of line 
flux, ionized gas velocity, and ionization state across the outflow. In Section |^, we describe our new 
data for the disk, halo, and outflow and provide comparisons with other published observations 
and a number of kinematic models. We present our conclusions in Section |5|. A subsequent paper 
describes refinements to our kinematic and ionization models to accommodate new observations 
from the Keck and Hubble Space Telescopes. 



2. Observations and Reductions 

2.1. Fabry-Perot Observations 

Observations of M82 were obtained in February, 1986 at the 3.6-meter Canada-France-Hawaii 
telescope (CFHT) on Mauna Kea, in the emission lines of Hq A6563 and [N II] A6583. The Hawaii 
Imaging Fabry-Perot Interferometer (HIFI; [Bland Sz Tully 1989 ) was used with an ET-50 etalon 



from Queensgate Instruments. The etalon coatings were formulated to provide a free spectral 
range of 85A and a finesse of 60, for a spectral resolution of 1.4A (65 km s~^, at Ha). The 
512 X 512 pixels of the Texas Instruments CCD, after 2x2 on-chip binning, each subtend 0'.'86 
(~13 pc at the distance of M82). This scale undersampled the estimated 1" seeing disk but was 
necessary due to the low flux levels of the outflow emission. The 3^5 field covers almost the entire 
region of minor-axis optical filaments. Other characteristics of the observing system are provided 
in Table along with their experimentally determined values. 

A second set of Fabry-Perot observations were obtained in March, 1992 at the University of 
Hawaii 88-inch telescope on Mauna Kea, in the emission line of [O III] A5007. Again the HIFI 
system was used, with an ET-70 etalon owned by Dr. Sylvain Veilleux. Although the free spectral 
range of this etalon (6lA) was smaller than that used for the Ha and [N II] observations, the 
finesse was essentially identical, as were the effective spatial and velocity resolutions (0'.'85 pix~^ 
and 60 km s^^ at [O III], respectively) and the field of view. The characteristics of the observing 
system are provided in more detail in Table ^. In both cases, the etalon was tilted to remove 
internal reflections and paired with wide (50A FWHM) interference filters to eliminate interorder 
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confusion problems. 



2.2. Data Reductions 



Most of the data reduction was performed under the Zodiac image processing system 
(Miyashiro 1982; Shopbell 1997), encompassing standard routines as well as a number of programs 
written by the authors and optimized for Fabry-Perot data. IRAF^ was also used for a small 



number of specific problems. A summary of the data reduction steps is now provided; see Shopbell 



1995 for more details. 



2.2.1. CCD Reductions 

The bias level was subtracted from each frame, using the average of several CCD bias frames. 
Cosmic rays were denoted by hand and removed by nearest-neighbor interpolation. Bad columns 
were also denoted by hand and then removed by sine interpolation across each row of the bad 
column. The sine function was selected for use in the interpolation because of its favorable Fourier 
properties: its Fourier transform is the top-hat function, implying even sampling of the image 
noise structure over a finite range of frequency. 

A flatfield image was constructed for our data by summing a series of frames of a white 
source, scanning the etalon through a free spectral range. Any first-order slope in the illumination 
pattern was removed with a two-dimensional linear fit. (Sky flats were not available.) The flatfield 
was then divided by a smoothed (o" 3 pixels) version of itself, normalized, and divided into each 
data and calibration frame. Due to problems with on-chip binning, the movement of transient 
defects between the images of each cube, and temperature fluctuations of the CCD, the flatfielding 
process did not substantially alter the noise content of most frames. Fortunately, the pixel-to-pixel 
sensitivity variations comprise a noise source at a level of only ~1.6%, a small value compared to 
other Fabry-Perot photometry uncertainties. 



2.2.2. Fabry-Perot Reductions 

Cube building. In order to randomize variations in the sky background due to changing hour 
angle, approaching dawn, etc., as well as instrumental fluctuations, the data frames were observed 
in a random or staggered order with respect to the etalon optical gap. The frames were first sorted 



^IRAF is distributed by the National Optical Astronomy Observatories, which is operated by the Association 
of Universities for Research in Astronomy, Inc. (AURA) under cooperative agreement with the National Science 
Foundation. 
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by gap and stacked to form cubes: target (M82) data cubes through each filter, a whitehght or 
flatfield cube through each filter, a standard lamp calibration cube, and standard star cubes. 

Whitelight calibration. The whitelight cube was used to map the spatial and spectral shape of 
the interference filter, which was then removed in the same manner as flatfielding. The whitelight 
cube was heavily smoothed {a ~ 15 pixels) in the spatial dimension, re-sampled from 0.65A to 
O.99A steps, to match the etalon gap spacings of the data cubes, then normalized and divided into 
the corresponding data cubes. 

Frame alignment. The frames in each data cube were aligned using two stars near the disk 



of the galaxy. (The single bright star in the field, AGK 3+69 428 [[Bettoni fc Galletta 198^ , 
2f5 southwest of the nucleus, saturated the detector.) Fractional pixel shifts were made by 
two-dimensional spline interpolation of the images. The spatial registration is accurate to better 
than one pixel Of! 5). 

We note here that, although the HIFI system and the tilted etalon provide a field that is 
relatively free of ghost reflections, the location of the bright star AGK 3+69 428 in the southwest 
corner of the frame was such that a pair of concentric ghost images of the star appear opposite 
the optical axis, in the southeast corner of the frame. Following frame registration, a region 
encompassing the ghost images was masked from further analysis. This is unfortunate, in that 
the minor axis emission of M82 runs directly through this region, but it demonstrates the care 
that must be taken when observing with Fabry-Perot systems, which always have prevalent ghost 
patterns. 

Ring fitting. Temperature and humidity variations produce drifts in the the expected spectral 
response of the Fabry-Perot etalon. In order to parametrize the spectral and spatial drifts in the 
Fabry-Perot system, we obtained a set of calibration lamp images, taken periodically throughout 
the night, all at the same etalon spacing. Elliptical fits to the rings revealed no noticeable 
variations in the radius or circularity of the rings, indicating that the spectral stability of the 
etalon was extremely good. Flexure of the telescope system was detectable as a two pixel (~1'.'5) 
shift of the ring centers (i.e., the optical axis) over the course of each night. The frame alignment 
procedure has removed this, with minimal effect along the wavelength axis. 

Data cube resampling and smoothing. While our Ha+[N II] data was sampled regularly 
at O.99A (45 km s^^) intervals across the Ha line, observing time constraints required us to 
interpolate several frames across the [N II] portion of the spectra, where the sampling was a factor 
of two coarser. Although a spline interpolation was performed without difficulty, a slightly larger 
error should be assumed for the final [N II] fluxes and velocities. The [O III] data set required 
extrapolation of a single (continuum) frame at either end of the spectra. 

A light Hanning filter was then applied along the spectral axis of each cube, in order to allow 
efficient automated fitting of the large numbers of the spectra. Tests indicate this had a negligible 
effect on the final fit parameters. 
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Phase calibration. The instrumental profile of a Fabry-Perot interferometer is a complex 
function of spatial position, wavelength, and optical gap spacing, given by the well-known Airy 
function ( Bland fc Tully 1989| ). Due to the large free spectral range and low interference order of 



the etalons, the monochromatic "phase surfaces," as observed in the emission lines of a calibration 
lamp, were well parametrized by the analytical expression for the three-dimensional Airy function. 
A fit to this function determined several system constants listed in Table |2|. This fit was then 
used to shift each spectrum in the data cube by the appropriate value to generate monochromatic 
frames. The convergence of two night sky lines into a single frame each confirmed the accuracy 
of the phase correction for the Hq;-|-[N II] data set. For the [O III] data set, however, a poorly 
sampled calibration cube prevented us from performing an accurate phase correction. Therefore 
the [O III] observations will be used for flux measurements and morphology, but not for kinematic 
studies. 

Sky subtraction. A limited field of view prevented us from obtaining a sky spectrum devoid of 
galaxy emission. We therefore removed the two bright night sky lines in the Ha+[N II] data set 
by subtracting Gaussian components with the proper velocity and a mean fiux level as observed 
across the field. The lines were identified as OH emission at 6553. 6lA and 6577. 28A ( Psterbrock 



Martel 19921) , providing the wavelength calibration for the spectral axis. The night sky continuum 



was not removed from the data. The [O III] data were not sky-subtracted, due to the low level of 
sky flux and the difficulties associated with the non-phase-corrected data. Resulting errors in the 
final spectral fits are negligible. 



2.2.3. Photometry 



Recent studies (e.g., Veilleux et al. 1994| ) have shown that Fabry-Perot data can be fiux 



calibrated to an accuracy equal to that of CCD imagery and longslit spectroscopy. Rather than 
scaling to an externally calibrated data set, we have applied a primary calibration, employing 
standard star observations directly. Since this Fabry-Perot fiux calibration procedure is described 
in more detail elsewhere (e.g., Bland-Hawthorn et al. 1997), we outline only the salient points 
below. 

The star e Orionis was used to calibrate the Ha+[N II] observations; a Lyrae and rj Hydrae 
were used for the [O III] observations. After cosmetic cleaning and applying the whitelight 
correction described above, the counts in each stellar image were summed, and then scaled by 
the telescope aperture size and exposure time. To obtain pixel values in units of counts cm~^ 
sec~^ we must then divide by the effective spectral bandpass of the etalon at each pixel. 
We emphasize that this value is not necessarily the same as either the spectral sampling of the 
cube or the etalon resolution. Rather, it is a measure of the amount of flux transmitted at a 
given wavelength for a sequence of etalon spacings; it depends primarily on the finesse {Nn) 
of the etalon. This "effective bandpass" was calculated by summing the flux under a synthetic 
monochromatic spectrum spanning an entire order of the appropriate theoretical Airy function 
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and dividing by the peak value of the Airy function. We derived an effective bandpass of 3.0A for 
our data sets. The scaled values of observed stellar counts were then compared with published flux 
values ( payes 1970| ; Hayes & Latham 1975), corrected for atmospheric extinction. Each pixel in 



the final spectra has units of ergs cm~^ sec~^ pixel"^ frame"^, yielding total fluxes from the line 
fits in units of ergs cm~^ sec~^ pixel"^. The final estimated systematic error in the flux calibration 
was ~7.5%. 



3. Fabry-Perot Maps 

The emission-line spectra were fit with Gaussian profiles across the field to generate spatial 
maps of line characteristics such as flux, velocity, and dispersion. Single-component fits were made 
to the Ha and [N II] lines across the majority of the field, well into the halo of the galaxy. There 
are regions north of the galaxy's disk where double components were fit to the Ha line profile, as 
well as regions south of the disk where double components were fit to both the Ha and [N II] lines. 
The [O III] line was fit with a single Gaussian component throughout the field, although visual 
inspection revealed signs of line splitting in specific regions. 

The region of acceptable fits is identical for the Ha and [N II] maps, although the [N II] fits 
in the northern- most regions of the maps have significantly higher errors, due to clipping of the 
6583A profile by the limited spectral coverage. The region with sufficient [O III] flux for profile 
fitting is much smaller than that of Ha or [N II] . The central regions of the galaxy are saturated 
in both the Ha and [N II] observations. 

We now describe the important features of each Fabry-Perot map. Each spatial map is 
presented at approximately the same scale, with tick marks separated by one arcminute (~950 pc). 
The maps from the Ha+[N II] data set contain regions in which the spectral lines are split; 
these maps show both a flux-weighted total for the galaxy as well as values for the individual 
components separately. These components are referred to as the high-velocity component (HVC) 
and the low- velocity component (LVC). North is up and east is to the left in all images. The 
position angle of the major axis of the galaxy is approximately 65°. 

Figure |^ illustrates the flux distribution in the light of Ha 6563 A and [N II] 6583 A. We 
measure a total unsaturated Ha flux of 9.9 x 10~^^ ergs cm^'^ sec~^ and a total [N II] flux of 
4.1 X 10~^^ ergs cm~^ sec~^. A rough comparison with the Ha imagery to be presented in the 
next section indicates that the saturated nuclear regions contribute an additional ~ 2.1 x 10~^^ 
ergs cm~^ sec~^ in the Ha line, or ~17% of the total. The flux is concentrated in the nucleus and 
along the minor axis, with very little emission originating in the extended disk of the galaxy. The 
nuclear line emission has saturated the detector in two concentrations ("the eyes"), as well as a 
number of weaker ancillary regions. In contrast, the minor axis emission is spatially extensive 
and filamentary. Numerous long radial filaments can be seen extending more than a kiloparsec 
from the nucleus, as well as complex small-scale structures. Note, for example, the bright bow 
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shock-like arc approximately 500 pc SSE of the nucleus. 

The morphology of the extraplanar gas differs between the two sides of the galaxy, appearing 
more chaotic on the north side and showing signs of collimation in the south. The southern Ha 
emission exhibits a sharp edge on the eastern side. Furthermore, the emission to the south is 
considerably brighter and more extensive, although there appears to be an abrupt reduction in 
flux at a distance of approximately 500 pc from the nucleus. Examination of the Ha and [N II] 
emission maps reveals that, beyond 500 pc, the flux in the HVC drops rapidly, while the flux in 
the LVC remains more uniform. We also find pervasive diffuse emission at a low level throughout 
the halo of the galaxy. 

Figure § is a map of the flux from M82 at [O III] 5007A. We measure a total [O III] flux 
of 2.6 X 10^^^ ergs cm~^ sec~^, more than an order of magnitude below that seen at Ha and 
[N II]. There is essentially no emission from the disk at this wavelength; the flux originates almost 
entirely within the nucleus and along the southern minor axis of the disk. Although the radial 
extent of the minor-axis emission is much smaller than that seen in [N II] and Ha, the filamentary 
morphology is similar. Indeed, excellent correlation is seen with the Ha flux map, which has been 
contoured on Figure |2[ A striking feature of the [O III] flux map is the presence of two distinct 
streams of emitting gas, each extending along the southern minor axis from one of the bright 
nuclear emission regions. These streams can also be identified in the Ha flux map (Figure ||), 
although the much more pervasive nature of the Ha emission makes the structure more difficult 
to discern at small radii. 

Figure |3| is a radial velocity map for the Ha-emitting gas in M82. Two prominent trends are 
seen: first, there exists a strong velocity gradient along the major axis of the galaxy, consistent 
with normal disk rotation, with the eastern end moving away from the observer. Second, a velocity 
transition is evident along the minor axis of the galaxy: strong blue-shifting of the Ha emission 
is seen south of the disk; strong red-shifting to the north. The major-axis rotation signature 
merges into this minor axis motion in a gradual fashion, although this may be due in part to the 
flux-weighting procedure. 

The HVC consists of highly redshifted emission in the north and highly blueshifted emission in 
the south {vr ~ 300 km s~^). The LVC consists of emission at roughly the systemic velocity of the 
galaxy {vgys ~ 203 km s~^; |de Vaucouleurs et al. 199ll ), in both the north and south. The northern 
and southern component pairs mirror each other almost exactly in terms of relative velocity and 
kinematic structure. Perhaps surprisingly, the velocities within each individual component exhibit 
little radial or azimuthal variation. The exception to this is the inner portion of the HVC, where 
the radial velocities approach the systemic velocity. 

A map of the line-of-sight velocities of the [N II] -emitting gas was also produced, but is not 
presented here, as it provides little additional information. The trends along both the major 
and minor axes are identical to those seen at Ha, but at a lower signal-to- noise ratio due to the 
decreased flux of [N II] line emission. A velocity map of the [O IIl]-emitting gas was not produced. 
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for the reasons outlined above. 

Figure |^ shows the logarithm of the ratio of the [N II] 6583A flux to the Ha 6563A flux from 
M82. The most remarkable feature in this map is the presence of a distinct region of low line ratio 
(0.0-1.0) along the minor axis, south of the disk. This region extends at least a kiloparsec in radius 
from the nucleus and is narrower than the regions of minor axis filaments visible in the Ha and 
[N II] flux maps (Fig. |^). The region of low ratios can be separated into two distinct structures, 
originating with the bright central eyes of M82, similar to the structures seen at smaller radii in 
the [O III] 5007A flux map (Fig. g). 

Elsewhere in the galaxy, higher line ratios 1.0) prevail, particularly at large radii within 
the inclined disk. Discrete H II regions can be seen in the disk as small concentrations with the 
expected [N II] /Ha ratio of ~0.5 QOsterbrock 1989|) . The regions of extremely low line ratio to the 



north are a result of the incompletely sampled [N II] line, and should be ignored. 

Figure |5| is a map of the logarithm of the ratio of the [O III] 5007A flux to the Ha 6563A flux 
from M82. The ratio has been calculated over the entire spatial extent of the [O III] emission, as 
seen in the flux map (Fig. Q). The observed ratios are low ('^0.05) throughout the nucleus and 
regions south of the disk, with the only apparent trend being a gradual increase in the ratio with 
distance from the nucleus. 



4. Discussion 

We now discuss the observational implications of these data in the context of each of the 
galaxy's primary morphological components: starburst disk, extended halo, and galactic-scale 
outflow. 



4.1. Starburst Disk 



Our observations reveal a small irregular disk containing heavy obscuration and a central 
starburst. The Ha and [O III] flux maps (Figs. |l] and ||) reveal little ionized gas in the extended 
disk of M82, outside of the central nuclear starburst region, as found by slit spectra ( p'Conneli 



fc Mangano 1978 ). Small concentrations of line emission that do exist in the disk outside of 
the starburst, such as the regions 0^5 east of the nucleus (see also p'Conneli fc Mangano 1978 



Fig. 3), appear to be giant H II regions. The [Nll]/Ha ratio map (Fig. ^) indicates moderate 
values for these regions (~0.56), consistent with a cooler range of photoionized H II regions 
(Tmax ~ 40,000 K, [Evans fc Dopita 1985| ; [Veilleux fc Osterbrock 1987| ). Such regions exhibit 
low levels of [O III] emission, consistent with our non-detection in Figure 0. Although the high 



extinction in the disk of M82 (3-27 mag; [O'Connell et al. 1995| ; [McLeod et al. 1993| ) serves to hide 
smaller H II regions from our optical observations, the paucity of large H II regions suggests that 
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the level of star formation must be low outside of the central starburst. This conclusion has also 



been reached by multi-epoch radio supernovae studies (e.g., Kronberg, Biermann, & Schwab 1985 



Huang et al. 1994). The weak nature of the [O III] emission implies a high overall metallicity in 



M82, as has been suggested by previous observations (e.g., Duffy et al. 1987| ; Paffney fc Lestei 
1992 ), and as is anticipated due to chemical enrichment by the starburst (e.g., Kobulnicky &j 



Skillman 19961) . 



The nuclear emission is dominated by two large saturated regions, each approximately 200 pc 
in diameter, and centered ~125 pc from the 2 ^m stellar nucleus (see Fig. ||). These regions 
correspond to knots A and C of O'Connell &: Mangano 197^ and are known to be extremely high 
surface brightness "clusters of clusters" of young (T ^ 50 Myr) stars (lO'Conneh et al. 1995^ . Also 
identifiable in Figure |l| are knots D and E (also saturated), as well as knots G, and H. Knot 
B is extremely faint at Hq wavelengths, especially when compared to broadband (p'Connell & 
Mangano 197^ ) and ultraviolet ( Hennessy 1996| ) observations, suggesting a higher gas content in 
the inner regions of the galaxy. The outflow can be traced to knots A and C. This relationship is 
particularly well-demonstrated by the [N II] /Ha map (Fig. Q) and the [O III] flux map (Fig. ^). 

We have plotted the Ha radial velocities along a straight line corresponding to the major axis 
of the galaxy in Figure ^ (panel c), along with a number of rotation curves from the literature 
at a range of wavelengths. The published systemic velocity of M82, 203 km s~^ ( |de Vaucouleurs 
et al. 199l| ), corrected to a heliocentric value of 208.7 km s~^, has been subtracted from the Ha 
Fabry-Perot data. The rotation curve rises to approximately 100 km s~^ within ~9" (~140 pc) of 
the nucleus and remains relatively flat to the edges of the observations. No substantial fall-off in 
the rotation curve is seen, due primarily to the limited radial extent of line emission from the disk. 
The observed rotation curve matches well those found at Ha by Heckathorn 197^ (his Fig. 12) and 



McKeith et al. 1993| (panel d of Fig. |6|), including the turn-over and asymptotic velocities and the 
nuclear velocity gradient (~11 km s~^ arcsec"^). Figure |6| also illustrates an increase in the central 



velocity gradient with wavelength. As has been pointed out by other authors (e.g., McKeith et al. 



1993), this effect is indicative of the large extinction toward the nuclear regions of the galaxy. 



The Ha rotation curve also correlates well with that of the 250 pc nuclear ring seen in 



molecular emission lines (e.g., Loiseau et al. 1990 ). The double-lobed structure of the central Ha 
emission is suggestive of an edge-on ring structure as well, interior to the molecular ring. The 
starburst could be identified with the ring of ionized gas and is probably propagating outward, 
fueled by the cold gas in the molecular ring ( [Waller, Gurweh, fc Tamura 1992| ; cf. |Shen fc Lc| 



1995 ). This ring may also provide much of the extinction toward the nucleus of M82 (e.g., Telesco 



fe Gezari 1992). The dynamic conditions of the central regions would have removed most of 



the obscuring material interior to the ring, as has been suggested in our Galaxy (e.g., Becklin 



Gatley, & Werner 1982). The resolution of the central star clusters by recent HST observations 
(O'Connell et al. 1995) suggests that any such ring would be very clumpy however, and the two 
bright regions may indeed represent real spatial enhancements in the distribution of ionized gas 
and star formation. The relationship between this ring structure and the proposed bar in M82 
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(r ~ 500 pc; fTelesco et al. 1991 ; |Achtermann &: Lacy 1995 ) is not clear. 

The [N II] /Ha line ratio map (Fig. ^) shows a high ratio in the disk, especially at large radii, 
where the value is observed to exceed 1.0. The actual line ratios in the disk may be even slightly 
higher and possibly more uniform, due to dilution by the central starburst and the halo. Such 
high ratios are by no means rare in the nuclei of disk galaxies ( [Keel 1983 ), particularly LINERS 
(e.g., NGC 4319 ^ulentic fc Arp 1987| , NGC 5194 [ [Ford et al. 198^ ; [Goad fc Gallagher 1985| ), and 
are understood to originate with shock excitation and/or photoionization by a power-law source 
( Veilleux fc Osterbrock 1987] ). Although it is unlikely that M82 harbors an AGN (e.g., Rieke 
et al. 19801; [Colina fc Perez-Olea 199^ ; [Muxlow et al. 1994[) , a comparison with other galaxies 



which exhibit pervasive high [N II] /Ha line ratios in their extended disks is instructive: studies of 
NGC 3079 ( [Veilleux et al. 1994[ ) and especially NGC 1068 ( [Bland-Hawthorn, Sokolowski, fc Cecil 



1991) have found high [N Il]/Ha ratios of 0.6-1.3 across much of the inner disk. In NGC 1068, 
the disk H II regions are found to reside in regions of lower [N II] /Ha ratio (~0.3-0.8), while the 
disk as a whole is permeated with gas exhibiting the higher ratios, similar to what we observe in 
M82. In order to boost this forbidden line to recombination line ratio, the "heating per ionization" 
must be high, for which most models require high energy photons, energetic electrons, or a dilute 
radiation bath ( Sokolowski 199^ ). 



In the case of M82, an attractive candidate to enhance the [N II] /Ha ratio in the disk is the 
concept of "mixing layers" ( Slavin fc Cox 1993| ; Voit, Donahue, &: Slavin 1994| ). The turbulence 



resulting from the interaction of hot supernovae remnants with the ambient ISM creates an 
intermediate temperature phase, which models suggest emits a radiation field somewhat harder 
than thermal bremsstrahlung. This process can produce [N Il]/Ha ratios in excess of 1.0 ( flavin. 



Shull, & Begelman 1993). By employing soft X-rays as the photoionization mechanism, mixing 



ratios have been used to model line emission ratios of up to 3-4 in cooling flows (Donahue &; 



Voit 1991; Crawford fc Fabian 199^ ). Considering the current optical appearance of the disk of 



M82, its interaction with the galaxy M81 approximately 10 years ago, and the energetic activity 
associated with the nuclear starburst, a turbulent inner disk would not be unexpected. 

However, as mentioned above, our observations do not detect significant numbers of star 
forming regions outside the starburst nucleus. Even if this is attributed to high levels of extinction 
in the disk, other wavebands confirm the low level of star formation. All of the radio supernovae 



discovered by [Kronberg, Biermann, fc Schwab 1985[ are within 300 pc of the galaxy's nucleus, well 
inside the region of highest [N II] /Ha ratios. The diffuse X-ray flux in the disk has also been shown 
to decrease rapidly with radius, implying a reduced star formation rate outside the nuclear regions 



(Bregman, Schulman, & Tomisaka 1995). Other models for producing high [N II] /Ha ratios in the 
disk must be considered, such as chemical enrichment and cosmic ray heating. A combination of 



shock and photoionization may also provide a solution (e.g.. Hunter &: Gallagher 1990 ), although 
detailed models are not yet available. 
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4.2. Extended Halo 



Early polarization observations detected a strong linear polarization throughout the halo of 



M82, with the position angles oriented perpendicular to a radial vector from the nucleus (Elvius 



1969 ). This polarization was interpreted as evidence for a scattering component in the galaxy (e.g., 



Solinger 19691) , probably consisting of electrons illuminated by the nucleus. The more difficult 
issue has been the polarization of the optical filaments themselves. The first observations in this 
regard ( Visvanathan &: Sandage 1972| ) determined that the minor axis filaments and the halo were 
equally polarized at optical wavelengths, suggesting that perhaps there was no "explosion" in M82 
and that the off-axis filaments were merely density enhancements in a dusty cloud through which 
the galaxy was moving ( polinger, Morrison, fc Markert 1977 ). 

Although the discovery of split emission lines ( Axon fc Taylor 1978| ) and the minor-axis X-ray 
halo ( Watson, Stanger, &: Griffiths 1984| ) have virtually eliminated this alternate interpretation 
for the optical filaments (although see [Rohan, Morrison, fc Sadun 198^ ), the polarization 



measurements remain poorly understood. Recent observations (e.g., Bcarrott, Eaton, & Axon 



1991 ) indicate that the optical filaments may indeed contain a scattered component, but even then 



there is uncertainty as to the source of the polarization, i.e. the nucleus ( Visvanathan 1974 ) or the 
entire disk ( ^olinger fc Markert 1975| ). Unfortunately, with few exceptions (e.g., Schmidt, Angel 



fc Cromwell T976[ [Bland fc Tully 1988| ; poane 1993| ; guchkov et al. 1994[ ), the presence of a halo 
component in M82 has been largely ignored. 

Our detailed analyses in Ha and [N II] clearly confirm the existence of the smooth exponential 
halo noted by Bland & Tully 1988. Figure |^ illustrates the flux along a narrow (~9") band parallel 
to and approximately 45" southeast of the major axis of the galaxy. The flux due to the outflow 
filaments can be seen superimposed upon an exponentially decreasing background. The halo has 
be detected across our entire region of fit lines, approaching a fiux level of 10~^^ ergs cm~^ sec~^ 
arcsec"^ at a radius of 1 kpc. 

We have estimated the radial profile of the halo flux at Ha with an exponential function: 

/(Ha) = /o(Ha)e-"/"% (1) 

where lo(Ha) ~ 1.6 x 10"^^ ergs cm~^ sec~^ arcsec"^ A^^ and re ~ 315 pc. This proflle has been 
overlayed on the cut in Figure ^. For the observed halo line width of -^350 km s~^, the integrated 
Ha fiux is ~ 2.5 x 10^^^ ergs cm~^ sec^^. This compares with an observed fiux from the filaments 
and unsaturated nucleus of ~ 9.9 x 10~^^ ergs cm~^ sec~^. For a distance of 3.25 Mpc this implies 
a total Ha luminosity from the halo of 3.2 x 10^^ ergs sec~^. 



The azimuthally symmetric polarization pattern of the halo in broadband light ( Schmidt 



Angel, fc Cromwell 1976| ) suggests a scattering origin. The halo is known to comprise cold neutral 
atoms ( [Cottrell 1977 ), relativistic electrons ( Seaquist fc Odegard 1991 ), dust ( [Visvanathan & 
Sandage 1972) and warm ions ( [Bland fc Tully 1988 ). If we associate the line-emitting halo with 
the polarized component, the line width (~350 km s~^ FWHM) reflects either the motion of 
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scattering mirrors embedded in a warm medium or the "beam-averaged" projected kinematics of 
the nuclear and large-scale disk gas. If the observed line dispersion arises from thermal motions 
of electrons, the kinetic temperature must be less than 1000 K, at which point the flux from 
recombination would overwhelm the scattered flux for any reasonable halo density {nn ~ 1 cm~^; 
Cottrell 1977| ). Dust scattering is expected to be much more efficient in any case. The ratio of the 



scattering optical depths can be written 



= (2) 



where na and rig are the dust and electron densities; ad and Ue are the dust and electron scattering 



cross sections. For ad, we form the more conservative weighted mean of the MRN (Mathis, Rumpl 



fc Nordsieck 1977 ) grain distribution for which a^, ~ 0.01/um. In the local ISM, the dust-to-gas 
number density ratio is nd/nu ~ 10~^^ ( Pstriker &: Silk 1973| ) although this depends strictly on 
the grain composition. With the conservative assumption of a totally ionized halo, we deduce a 



scattering ratio Rde ~ 50, verifying the relative importance of dust scattering. Schmidt, Angel, &: 



Cromwell 1976 demonstrate that a reasonable optical depth to dust scattering is ~ 0.1 from a 



comparison of the halo broadband flux and the estimated disk flux. 

The halo dust could reside in an extended neutral or warm ionized medium, some fraction 
of which could be supplied by the energetic wind ( [Burbidge, Burbidge, fc Rubin 1964 ). We now 



compare the timescale for dust destruction by sputtering with the estimated age of the starburst 
wind, Tyjind ~ 3 X 10^ years (e.g., Lynds &: Sandage 1963 ; Bland fc Tully 19*8^ ). From Ostriker fc 



Silk 1973, the timescale for grain sputtering is 

Tsput ~ 3 X 10^ ne (^^) years, (3) 

where Y is the sputtering "yield," i.e., the number of atoms released per impact. Y varies from 
~ 5 X 10^"^ at the threshold sputtering temperature of 3 x 10^ K to 0.01 for temperatures of 
10^-10^ K. 

At the low temperatures expected in a galactic halo, T <^ 10^ K, the lifetime of dust is 
greater than 10'' yrs, comparable to or longer than the lifetime of the outflow. However, any dust 
located directly in the hot (T ~ 10* K) X-ray-emitting wind survives for no more than 10^ years, 
and has therefore been destroyed. Such an effect is supported by low-resolution radio maps of 
the M81/M82 region, which indicate an anti-correlation between the wind lobes and H I column 
density ( [Cottrell 197^ ). We note that more recent studies of gas-grain sputtering and grain-grain 
collisions (|Tielens et al. 1994 ; Jones, Tielens, &: HoUenbach 1996 ) suggest that grains may be able 



to survive much longer than previously thought. However, most of these studies are specific to 
the three-phase ISM in the Galaxy, and it remains unclear precisely how the results should be 
extended to galactic wind systems, which characteristically involve higher temperatures (10*^ K vs. 
10^ K) and larger velocities (600 km s~^ vs. 200 km s~^) than standard ISM models. 
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Given the sputtering timescales, if the dust has been deUvered into the halo by the wind, 
it must be as a component of cooler material entrained by the hot wind itself. However, this 
implies that the halo dust and optical emission line filaments may then have the same origin and 
morphology, yet the polarization observations do not seem to indicate a minor-axis concentration 
in the dust distribution ( ^olinger fc Markert 1975 ; Schmidt, Angel, & Cromwell 1976; Scarrott, 



Eaton, fc Axon 1991 ). We also do not observe substantial redshifted emission south of the galaxy, 



as would be expected from mirrors moving with the outflow. 

Alternatively, dust may have been forced into the halo at an early stage of the outflow, at 
a time when it was dominated much more by radiation from massive stars in the central burst 



than by supernovae. Such a mechanism has been hypothesized by Ferrara 1997 to explain the 
appearance of high-z dust in the Galaxy and other edge-on spirals. It has been demonstrated that 
radiation pressure is sufficient to evacuate a large fraction of the dust near an active star-forming 
region into the halo, creating a dust distribution which varies slowly with height above the disk. 
Given the current importance of radiation effects in the inner portion of the M82 outflow (see 
Section 4.3.6| below) and the extensive nature of the observed dusty halo, such a scenario seems a 
reasonable model. 

Regardless of any mechanism of relocating disk dust into the halo, clearly the hypothesized 
encounter between M81 and M82 ~ 10^ years ago (|CottreU 19771 ; |Yun, Ho, fc Lo 1994|) has played 



an important role in the evolution of the halo in M82. (This encounter is also thought to have 
initiated the central starburst in M82, e.g., Mihos &: Hernquist 1994 ) Radio observations have 
shown massive clouds of H I surrounding both galaxies, with large arcs and bridges joining them 
and the nearby galaxy NGC 3077 (e.g., pavies 1974| ; |Yun, Ho, fc Lo 1994| ). The large-scale 
velocity structure of this gas blends with the global Ha velocity trends in M82, matching the 
systemic velocity and even the minor-axis velocity gradient ( Pottrell 1977] ) . The H I cloud is 
clearly extensive enough to replenish dust that has been destroyed by sputtering. Moreover, this 
massive reservoir of atomic gas should help to maintain the tenuous halo gas itself in M82, a 
component that is required by hydrodynamical models in order for the outflowing wind to produce 
observable structures ( Suchkov et al. 19941) . 



4.3. Galactic-scale Outflow 



A number of studies have concentrated on the minor axis filaments in M82. Imaging 



observations (e.g., [Lynds fc Sandage 1963| ; [Ichikawa et al. 19941) have been directed toward 



understanding the morphology of the outflow, while spectral studies (e.g., [Burbidge, Burbidge, fc 
Rubin 19641 ; peckathorn 1972| ; [Axon fc Taylor 1978| ; [Wilhams, Caldwell, &: Schommer 198^ ; [Bland 



fc Tully 198£ ; McKeith et al. 1995 ) have attempted to parametrize the kinematics of the filaments. 
The most recent observations detect split emission lines that suggest a pair of expanding bubbles 
or cones. The exact size and shape of the expanding structures remains a topic of some debate. 
Additional optical spectroscopy has measured the radial variations of line emission along the 
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outflow (e.g., Heckman, Armus, fc Miley 1990| ). Ratios such as [N II] /Ha increase with distance 
from the nuclear starburst, suggesting a gradual change in the gas excitation mechanism. We now 
investigate these findings in light of the spatial and spectral coverage provided by our Fabry-Perot 
spectrophotometry. 



4-3.1. Optical Morphology 

Optical emission line images of the M82 outflow reveal a complex distribution of radial 
filaments and knots along the minor axis of the galaxy. In order to examine the morphology of 
these filaments, R. B. TuUy kindly obtained for us very deep {t ~ 3000 sec) Ha imagery at the 
University of Hawaii 88- inch telescope on Mauna Kea (Fig. |8|). We imaged a 6' circular field, a 
factor of two larger than our Fabry-Perot field. The optical filaments associated with the minor 
axis outflow are visible across the entire field, to radial distances of at least 3.5 kpc to the north 
and 2 kpc to the south. This represents approximately half the extent of the soft X-ray halo (see 
below). The fllaments are bright within a kiloparsec of the disk, particularly in the south. In 
these inner regions, a web of long (~300 pc) filaments appear distributed roughly parallel to the 
minor axis, amidst a "foam" of smaller emission structures. Beyond a kiloparsec in radius, the 
filamentary structure breaks up into fainter, more isolated clumps. 

The faintness of the northern filaments is almost certainly due to obscuration by the inclined 
disk of the galaxy (cf. [Hennessy 1996 ). The inclination ang le of M82 is estimated to be 81? 5 (^ynd| 



fc Sandage 1963 ), such that the nearer edge of the galaxy is projected on the northwest side of 
the nucleus. We observe the nuclear regions of the galaxy through the southern side of the disk, 
and therefore expect the inner filaments to be much brighter there. If we assume a line-of-sight 



dimension for the optical disk equal to its linear dimension on the sky (~llf2; |de Vaucouleurs ei 



al. 199 1| ), we expect the northern outfiow to be at least partially obscured to projected radii of 



~ 0f8. This is consistent with the similar morphologies of the filaments beyond this radius in the 
north and south in Figure ^. 

Although complicated by obscuration in the north, the inner filament structure also appears 
to differ between the two lobes in terms of collimation. The bright inner kiloparsec of the southern 
lobe is relatively well confined to the minor axis, whereas the northern outflow fllaments cover a 
much wider range of opening angle. This is seen particularly well in the Ha flux map from the 
Fabry-Perot data (Fig. |^), and is probably not an obscuration effect, but rather a difference in the 
physical morphology of these two lobes. One possible interpretation is that the nuclear starburst is 
located slightly above the galactic mid-plane. The smaller mass of covering material to the north 
would make collimation of an expanding wind more difficult, resulting in an almost immediate 
"breakout" of the wind from the disk. 



Finally, we have noted that the [N II] /Ha ratio map and [O III] flux maps indicate that the 
southern outflow involves two distinct components, each originating from one of the central bright 
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emission regions. This may be a "limb brightening" effect as has been suggested for the pair 
of central emission regions themselves. This implies that the emitting filaments are distributed 
along the outer surface of the outflow, rather than throughout the volume. On the other hand, it 
is likely that the two bright central regions from which the outflow streams appear to originate 
are stellar "superclusters" (O'Connell et al. 1995) which merely happen to be physically located 



on either side of the kinematic center of the galaxy, from our point of view. Regardless, we do 
not observe emission enhancement along the outflow axis, as would be the case for a volume 
brightened distribution. 

If we do not include emission arising within approximately 8" (125 pc) of the disk, our line fits 
encompass an Ha flux of 7.6 x 10^^^ ergs s^^ cm^^. After a rough subtraction of the halo model 
given in the previous section, this implies a total filament luminosity L{Ha) ~ 7.6 x 10'*° ergs 
s~^ (cf. Tab. H). Assuming the filaments are completely ionized, using a case B recombination 
coefficient for T ~ 10,000 K of = 2.59 x 10"^'^ cm^ s~^ ( psterbrock 1989| ), and employing 



the outflow geometry to be discussed in Section 4.3.4 below, we calculate an rms filament density 



(ng) ~ 5 • /^^/^ cm~'^. A very rough estimate for the filament filling factor of / ~ 0.1 would suggest 
a mean electron density in the filaments (ng) ~ 15 cm~^, easily consistent with [S II] doublet 
ratios in the low-density limit throughout much of the large volume of outer filaments. This mean 
density implies a filament mass M ~ 5.8 x 10^ Mq distributed in a volume ~ 1.1 x 10® pc^, and 
a kinetic energy in the ionized filaments KE ~ 2.1 x 10^^ ergs. These values are consistent with 



those first computed by [Lynds fc Sandage 1963| , 5.8 x 10*^ Mq and 2.4 x 10^^ ergs, respectively. 



The entire mass of outflowing gas is estimated to be a couple orders of magnitude larger (Heckman 
Armus, fc Miley 199C| ). Note that the kinetic energy in the filaments is only about one percent of 
the estimated input supernovae energy (~ 2 x 10^^ ergs; Watson, Stanger, h Griffiths 1984). 



4.3.2. X-Ray Morphology 



A large x-ray halo has been observed in M82 (Watson, Stanger, & Griffiths 1984; Schaai 
et al. 1989t [Tsuru et al. 199(\ ; pre gman, Schulman, &: Tomisaka 1995| ), extending at least 5-6 



kiloparsecs in radius along the minor axis of the galaxy. Originally believed to be thermal 
emission from the hot (T ~ 10''-10® K) gas within the outflow, current models suggest an origin in 



wind-shocked clouds. J.N. Bregman has kindly provided us with a recently published (Bregman 



Schulman, &: Tomisaka 1995 ) image of M82 taken with the Rontgensatellit {ROSAT; piViimper 



1984D High-Resolution Imager (HRI; [Pfeffermann et al. 1987|) . The pair of 13' square images with 
total exposure time 33.8 kiloseconds were taken in early 1991 and late 1992. The resolution is 
about 5". 

Contours representing the ROSAT soft x-ray flux have been overlayed on our deep Ha image 
(Fig. P). The morphologies are markedly similar, with specific filaments and knots clearly enhanced 
in both emission bands. (Note, for example, the two extensive filaments in the north, and the 
knots approximately 500 pc SW of the nucleus.) But in light of the more extensive distribution 
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of the X-rays versus the optical emission, it seems unhkely that all of the X-ray flux is associated 
with specific filaments of optical emission. Since the regions of high optical/X-ray correlation 
are also some of the brightest in Ha, these knots presumably represent density enhancements in 
the outflow bubbles. The large extent of the X-ray halo and close spatial correlation with the 
Ha emission suggest an interpretation in terms of shocks driven by a fast, rarefied wind plowing 
into denser halo gas (e.g. Sutherland, Bicknell, fc Dopita 199^ ; Dopita &: Sutherland 1996 ). Our 
observed Ha/X-ray luminosity ratio of ~30 in the brightest filaments, identical to that derived 
by |Perez-01ea fc ColinaT996 , suggests a similar scenario. We defer detailed modeling of the Ha 
and soft x-ray emission to a subsequent paper, in which the shocks will be constrained by deep 
spectroscopy and imaging from the Keck and Hubble Space Telescopes. 



4-3.3. Optical Kinematics 

Our Fabry- Perot observations provide velocity information across the entire outfiow in M82. 
Two well-known features of this emission line gas should be emphasized a priori. First, the 
velocity signature of outfiow is clearly observed along the minor axis in both Ha and [N II] (see 
Fig. ^). The well-established disk inclination angle of 81? 5 ( Lynds fc Bandage 1963| ) and the 



increased opacity toward the northern filaments clearly establishes the large-scale kinematics as 
arising from outfiow (cf. [Elvius 1972| ; polinger, Morrison, &: Markert 1977| ). Second, the optical 
emission lines split into two distinct velocity components in large regions along the minor axis of 
the galaxy. The velocity structure of the individual components, as well as the difference between 
them, provides convincing evidence for the presence of outfiowing gas along the minor axis of 
M82 ( Axon &: Taylor 1978 ). Figure |9| illustrates the velocity structure of the Ha-emitting gas in 



M82, as observed in the frames of the Fabry-Perot data cube. Both minor-axis outfiow lobes span 
velocity ranges exceeding Vproj ~ 300 km s~^, overlapping at the systemic velocity. 

We now describe our attempts to understand the complex dynamics of the outfiow in M82 
with the aid of kinematic models. We have extracted one-dimensional velocity cuts and synthetic 
two-dimensional spectra from the Fabry-Perot data cubes. These cuts have been made relative 
to the outfiow axis, at an estimated position angle of ~150° (~85° relative to the major axis of 
the galaxy; [McKeith et al. 19951 ). While some authors have argued for a spherical wind in M82 



(e.g., [Seaquist fc Odegard 1991 ), most recent models find evidence for a bipolar outfiow, usually 
in the shape of cones with opening angles much less than 90°. This is borne out by our study: 
the line fiux, splitting, and ratio maps all exhibit marked azimuthal variations indicative of an 
aspherical outfiow morphology, strongly weighted toward the minor axis of the galaxy. Studies of 
the supernova distribution (e.g., [Kronberg, Biermann, Sz Schwab 1985 ) and other disk observations 



(e.g., Achtermann Sz Lacy 199^ ) indicate that even the central injection zone is probably not 



spherical, but rather in the form of a fiattened disk, with dimensions ~600 pc wide and ~200 pc 
thick. 



The specific velocities predicted by phevalier fc Clegg 1985 



do not agree with our observations. 
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An estimated supernova rate of 0.3 yr~^ ( |Rieke et al. 1980 ) yields an outflow velocity from their 
model of 2000-3000 km s~^ at the edge of the starburst injection zone. This value is several 
times larger than the deprojected velocities observed at optical wavelengths {v ~525~655 km s~^; 



see below) on much larger scales. However, the extensive radio continuum halo (Seaquist, Bell 



fc Bignell 1985| ; peaquist fc Odegard 1991] ) could arise from synchrotron radiation due to a 
population of relativistic electrons, as they are transported outward in a wind at velocities in the 
range considered by Chevalier & Clegg. The relationship between such a wind and the slower, 
denser minor-axis outflow seen at optical and x-ray wavelengths remains unclear. We suspect 
that much of the Ha filamentation arises from large-scale shocks from a high-speed wind plowing 
into the gaseous halo and entrained disk gas. Only a small fraction of the total wind energy is 



encompassed by the radio halo (~2%; Seaquist &: Odegard 1991 ) 



In modeling the bipolar outflow in M82, we model the individual line components rather than 
flux- weighted velocity profiles. In Figure we show the velocities of line fits to the dual Ha 
components along the minor axis of the galaxy. While the northern outflow components show a 
relatively constant projected separation of ~300 km s~^, the southern region of split lines reveals 
an intriguing variation. Within 200 pc south of the nucleus, where the favorable inclination of the 
disk allows us to measure line profiles closer to the starburst, the individual components of Ha are 
separated by a much smaller velocity, comparable to our resolution (~50 km s~^). Between 200 
and 500 pc from the nucleus, the components rapidly diverge, remaining at a constant separation 
of ~300 km s~^ beyond 500 pc. Maps of the line component splitting reveal a separation of 
this order throughout the spatial extent of both lobes. Also plotted in Figure 10 are fits to the 
Ha profiles from the long-slit optical spectra of McKeith et al. 1995| , which clearly confirms our 
observed minor-axis trend in the line splitting. An identical radial trend is seen in the [N II] line 
components along the southern outflow axis, but is not shown in the figure. 



4.3.4- Kinematic Models 

In order to understand the intrinsic velocity structure underlying the observed kinematics, we 
have undertaken two separate sets of three-dimensional Monte-Carlo wind simulations: a rounded 
expanding bubble and a pair of cones arranged as a funnel. Although primarily geometric in 
nature, these models allow us to estimate the three-dimensional morphology of the outflow, as well 
as intrinsic gas velocities, both of which are important constraints for physical models of the wind 
emission mechanisms and other observational phenomena. 

For our initial model, we used rounded bubble geometries, given by the spherical functions, 

p = A cos m9 — < 9 < +— (4) 

2m 2m 

p = A co^re -| <^<+| (5) 

(see Fig. p^g). Both functions produce a parabolic leading surface becoming conical at small 
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radii. The parameter m can be related to an opening angle and a Mach number in cases where 
the leading surface is produced in a bow shock. Truncating the bubble at inner and outer radii 
allows for the selection of a outflow structure with specific opening angle and radial curvature. 
We superimposed two alternate velocity laws upon this spatial geometry: a radial velocity vector, 
which corresponds physically to the case in which each gas parcel is accelerated by the flow 
originating at the bubble apex, or a velocity vector perpendicular to the surface of the bubble, 
which may be more appropriate for an expansion due to increasing heat and pressure inside the 
bubble, such as an inflating balloon. 

We then executed a series of Monte-Carlo simulations for each velocity law, varying bubble 
parameters such as opening angle, inclination angle, and inner and outer truncation radii. 
The conclusion reached from this set of simulations is that the rapid divergence of the velocity 
components in the southern outflow cannot be reproduced by a single bubble, at least not without 
invoking highly contrived velocity profiles for the wind. Studies of the minor axis x-ray distribution 
are similarly unable to model the emission with a single bubble or cone (e.g., puchkov et al. 1996 ). 



Morphologically, one can divide the southern region of split lines into two separate velocity 
regimes: the region within approximately 200 pc of the nucleus, where the split line components 
are separated by ~50 km s^^, and the region beyond 500 pc radius, in which the components are 
separated by a much larger, but still relatively constant, projected value of ~300 km s~^. The 
inner component is not observed north of the galaxy, presumably because the split lines cannot 
be resolved at a sufficiently small radius due to the intervening inclined disk of the galaxy. Note 
that the Ha and [N II] flux maps (Fig. |l|) also suggest the presence of two distinct regions, as 
the line flux drops sharply at the same radius at which the velocity components rapidly separate 
(~500 pc). 

We therefore performed another set of Monte-Carlo simulations, this time using a double cone 
geometry, given by the cylindrical function 

ni{z-ZQi), < z < 350 pc 

«^2(z - 202), 350 < z < 800 pc ^ ' 

(see Fig. |ll|6). The k parameters determine the opening angles of the cones, while the zq 
parameters control the radial extent of the cones through truncation. Again, we superimposed 
both radial and normal (i.e., perpendicular) velocity laws, but determined in the end that a 
velocity vector tangential to the cone surfaces provided the best match to the observations and 
was most easily understood physically, in terms of material entrained by the high-velocity wind. 
An additional parameter was used to smooth the abrupt projected velocity transition where the 
two cones meet (at z = 350 pc). 

After performing simulations over a range of cone sizes, inclination angles, opening angles. 



and velocity laws, we derived the final model shown in Figure 12. Parameters of the cones which 
best fit the observations are given in Table ^. The inner "cone" is almost a cylinder of radius 
equal to the injection zone (r ~ 200 pc), with a small, but non-zero, opening angle. The outer 
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cone has an opening angle of approximately 25°, in relative agreement with the "small" opening 
angle models of [Bland fc Tully 1988| and [McKeith et al. 1995 . Models in the "large" opening 



angle regime (e.g., 60°; Heckman, Armus, fc Miley 19901 ) do not match the observations, requiring 



excessively low intrinsic velocities and a larger projected spatial extent for the outflow region 
(see Fig. |l^). A large opening angle for the primary outflow cone also produces substantially 
skewed Doppler ellipses in spectra perpendicular to the outflow axis. This effect is due to the slit 
cutting through the back and front of the cone at different nuclear radii, and is not observed in 
our synthetic spectra (see Figs, and 12 ri). 



The inner and outer cones are inclined toward the observer by ~5° and ~15°, respectively. 



roughly aligning the back sides of the two cones (see Fig. lib). This is required to explain the 
lack of a sharp velocity gradient in the low-velocity component (LVC) at a radius of 350 pc, as is 
observed in the high- velocity component (HVC). In addition, since the LVC exhibits small but 
non-zero projected velocities, the back sides of the cones must be at a slight angle to the plane 



of the sky. These inclination angles agree with previous estimates (e.g., Burbidge, Burbidge, & 



Rubin 1964 ; Hennessy 1996| ). While the observed velocities of the HVC could be duplicated with 



smaller cone opening angles and a larger inclination, the low velocities of the LVC require small 
inclination angles. 

Initial attempts to model the kinematics of the outflow with a constant velocity law, i.e., 
using only the simple double-cone geometry to reproduce the observed velocity structure, were 
not successful. Figure ^ illustrates that both velocity components in the south and north exhibit 
non-zero slopes in the position-velocity plot. This can be understood as either a continuous 
change in the intrinsic gas velocity or as a change in the projected velocity through a continuous 
change in the outflow cone geometry. The latter case implies that both sides of the outer cone 
are constantly bending toward the observer, producing a slowly increasing projected velocity with 
radius. A gradually increasing intrinsic wind velocity is probably necessary as well, and can be 
understood from a physical standpoint. Buoyancy effects in the hot wind from the decreasing 
disk density with scale height, decreasing wind densities from the lack of collimation at larger 
radii, and other effects contribute to produce wind velocities that increase with radius in standard 
galactic wind models (e.g.. Chevalier fc Clegg 1985 ; ^uchkov et al. 1994| ). After testing a number 



of stronger power-law expressions for the gas velocity dependence on radius, we finally chose the 
simple linear model given in Table ^. Together with a constant cone opening angle, this intrinsic 
velocity structure produces linear projected velocity gradients that correspond well to those 
observed in Figure [l^. The intrinsic velocities of the gas range from 525 km s~^ near the nucleus 
to 655 km s^^ at a radius of 1 kpc. These velocities are comparable to the escape velocity for M82 
(see Table |l|), implying that the most distant entrained filaments are not bound to the galaxy. 
This conclusion is also supported by the large radial extent of the fast wind itself, as seen in soft 
X-rays, along the minor axis. 



Just as the outflow cones are not aligned with the minor axis of the galaxy along the line of 
sight, neither are they aligned in the plane of the sky. The region of split Ha lines constitutes 
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a cone on the sky with the expected opening angle of ~25°, but with a position angle of ~165°, 



approximately 15° greater than what the literature (e.g., McKeith et al. 1995) had previously 
defined as the "outflow axis." The cone axis is rotated ~100° from the major axis of the galaxy, 
and places the eastern edge of the cone almost directly parallel to the minor axis. In fact. 
Figure illustrates that this eastern edge is quite pronounced, both in the split emission lines 
and in the Ha flux observed from the inner collimated zone. This suggests that the tilting of the 
outflow cones in the plane of the sky has been produced by a relative density enhancement in the 
eastern lower halo which maintains collimation of the wind even as it fans out toward the west 
and toward the observer. 

While the large outflow cones appear to originate east of the galaxy's minor axis, a small 
region of split lines is also observed on the western edge of the collimated zone, approximately 
300 pc from the nucleus. The morphology of the Ha and [O 111] flux maps (Figs. || and |2|) indicates 
that this region constitutes a small bubble on the side of the larger outflow structure. We clearly 
see an enhanced rim of Ha emission around the bubble, and split lines within it. 

Recalling again the identification of two outflow "streams" in the [O 111] flux map and the 
[N ll]/Ha ratio map (Figs. ^ and|^, one might expect a more substantial outfiow from the western 
half of the nucleus, or at least a more centrally-positioned outflow cone. However, it appears that 
gas densities on the western side of the lower halo are substantial enough to keep the western 
stream from expanding into a cone. The stream appears to bend toward the west in the [O 111] 
flux map (Fig. |2|), and the only expanding structure that we observe is the one small bubble. 

But at larger radii, the ambient density toward the west must drop relative to the eastern 
side. While the outfiow remains tightly confined to the east, even beyond the collimated zone, 
certain Fabry-Perot maps show evidence of variations in the western side of the outer outfiow cone, 
suggesting a more azimuthally-extended morphology there, e.g., break-out into a less dense region. 
For example, the Ha velocity map (Fig. |3|) reveals high velocity clumps at radii of ~1 kpc that are 
distributed westward from the sharp eastern edge over almost 90° in azimuth. These clumps are 
clearly associated with the outfiow, and can also be seen as the blue-shifted emission just west of 



the split line region in our most distant outfiow cut. Figure 12d. The fact that this gas exhibits 
smaller velocities in panel c of this figure illustrates that kinematical effects from disk rotation are 
greater closer to the disk. In contrast, the emission immediately east of the split line region lies at 
effectively the same velocity in both panels c and d of Figure |l^, as this coincides with the sharp 
collimating edge on that side of the outfiow. The spatial structure of the Ha velocity map (like 



that published in Heckathorn 1972 ) also suggests the presence of rotating disk material that has 



been entrained and gradually diverted to the outfiow. 

To summarize, the inner 350 pc of the outfiow constitutes a fiow down a pipe. The outfiow is 
collimated, presumably by ambient and entrained disk material, and highly inclined to our line of 
sight, such that the observed radial component of the fiow velocity is only ~50 km s~^. Beyond 
350 pc, however, the collimation weakens and the outfiow expands rapidly as a cone of emission 
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with an opening angle of 25° and a projected front-to-back velocity separation of approximately 
300 km s-^ This expansion is preferentially toward the west and toward the observer. A linearly 
increasing intrinsic gas velocity with an initial value of 525 km and a gradient of 0.13 km s~^ 
pc~^ matches the observations well out to a radius of a kiloparsec. 

The line-splitting phenomenon indicates that the Ha-emitting filaments are produced on the 
surface of the outflow cones, at the interface between the wind and the ambient halo material. 
In addition to the minor axis velocity structure, this model explains the increased flux within 
the collimated region (see Fig. as a result of an elevated density with respect to the outer 
expanding cone {J- oc Tig). An increased density in the innermost regions has also been indicated 



observationally by the "filling in" of the line profiles at those radii (e.g., McKeith et al. 1995 ) 



While the inner bubble has an extremely small opening angle, we find that an outer cone of 
opening angle ~20°-30° fits the data most closely. This range agrees well with the "small" opening 
angle values found by other authors (e.g., ~30°; McKeith et al. 1995| ). Our observations do not 



support the "large" opening angle regime (e.g., ~60°; e.g., Heckman, Armus, &: Miley 199C| ). 



4-3.5. Comparison of Kinematic with Hydrodynamic Models 



Striking support for the two-zone model is obtained from a comparison with recent 



hydrodynamic outflow models. Initial two-dimensional simulations (e.g., Tomisaka & Ikeuchi 



1988) proposed an evolution of the outflow which included collimation of the flow at early stages 
by disk material, subsequent breakout of the confined flow along the minor axis, and eventual 
extension of the hot wind material up to a few kpc from the nucleus. 



More recent simulations ( Tomisaka &: Bregman 1993 ; Suchkov et al. 1994 ; Suchkov et al. 1996 ) 
include a separate corotating halo component and are able to reproduce much of the observed 
morphology through the interactions of the wind with the disk and halo. These simulations show 
that the outflow entrains disk gas around itself, dragging the cooler, denser material up to a 



couple of kiloparsecs above the plane of the galaxy (e.g., see Fig. 6 in Suchkov et al. 1994 ). The 
regions of densest entrained disk material, near the base of the outflow, serve to collimate the 
outflow beyond the height of the disk itself. The scale height of this collimation in the simulations 
is similar to that seen in our observations, ~500 pc. The "fingers" of disk material entrained 
to heights above the collimated zone can be identified with the optical emission line filaments 
we observe on the outer walls of the outflow cone. This entrained gas has also been observed at 
molecular wavelengths (e.g., Stark & Carlson 1984; Nakai et al. 1987; Sofue et al. 199d| ). Within 
the collimated region, the simulations show that both the wind and the confining walls are at 
their densest, consistent with the observed increased levels of optical emission. Finally, it should 
be noted that a recent analysis of the minor-axis X-ray distribution implies a partially confined 
outflow of this gas as well, within 1.6 kpc of the disk ( [Bregman, Schulman, fc Tomisaka 1995| ). 
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4-3.6. Comparison of Optical and X-ray Observations 
A number of authors have observed an extensive nebula of soft X-ray emission along the 



minor axis of M82 (e.g., 


Watson, Stanger 


& Griffiths 1984 


; Schaaf et al. 1989; 


Tsuru et al. 199C; 


Bregman, Schulman, & Tomisaka 1995|). 


The initial interpretation was that these are thermal 



X-rays, produced by the hot (T ~ 10^ K) gas in the outflow. The optical line emission must then 
arise from the cooler boundary layer of the wind, where the hot gas interacts with entrained disk 
and ambient halo material. Such a hot gas would not be gravitationally bound to the galaxy 



and could easily expand to the the observed radial distances of 5-6 kpc (Bregman, Schulman, fc 



Tomisaka 1995| ) in the estimated age of the starburst {t ~ 5 X 10^ yr; |Rieke et al. 1980| ; ^ oane &| 



Mathews 1993| ). Recent ASCA observations ( p?suru et al. 1994 ; Moran fc Lehnert 1997| ) find three 



temperature components in this wind, with more extended softer components, confirming that the 
hot wind is cooling as it expands. 

Recent hydrodynamic simulations and studies of the x-ray halo spatial distribution have cast 
doubt on this interpretation, however. Current hydrodynamic simulations ( [Tomisaka fc Bregman] 



1993 ; Buchkov et al. 1994 ) derive gas temperatures an order of magnitude larger than earlier 



estimates (e.g., Watson, Stanger, fc Griffiths 198'| ). This 10* K gas does not emit as strongly 



in the X-ray bands, producing insufficient soft and hard X-rays to account for the observations. 
Even Chevalier fc Clegg 1985| admit that the density falls too rapidly in their simple model to 



account for the x-ray photons as thermal. X-ray spectral studies derive a range of thermal gas 
temperatures and suggest alternative emission mechanisms ( [Fabbiano 1988| ; ^chaaf et al. 1989 ; 



Strickland, Ponman, fc Stevens 1996| ). Finally, correlations between our deep Ha imaging and 



Fabry-Perot observations and high-resolution ROSAT imagery support a non-thermal origin for 
at least a portion of the minor-axis x-ray emission. 



Figure 13 compares the spatial distribution of soft X-rays observed by ROSAT with our 
Fabry-Perot Ha flux map. As was evident at larger radii in the comparison with the deep Ha 
imagery (Fig. |^, the X-rays and optical line emission are clearly correlated. On large scales, the 
minor-axis x-ray flux drops at a radius of approximately 500 pc from the nucleus, as does the Ha 
emission. But also on scales as small as 10 pixels (150 pc), the x-ray and optical emission appears 
well correlated. This implies that soft X-rays are being produced in regions very close to those 
which are producing Ha emission, a situation which is very difficult to understand in terms of a 
thermal emission mechanism. 

These observations lend support to the x-ray emission mechanism suggested by several 



authors (e.g., Chevalier fc Clegg 1985; Suchkov et al. 1994; Strickland, Ponman, fc Stevens 1996): 



the soft X-rays arise from shocked disk and halo "clouds." This shocked gas can produce both 
the observed x-ray and optical emission, accounting for the spatial correlation in Figure |l3|. A 
hybrid model seems necessary in which the higher gas temperatures and densities near the nucleus 
create a region dominated by shocks at interfaces with disk and halo gas clouds, while the cooler 
temperatures and lower densities at larger radii produce a decrease in optical emission and an 
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increase in thermally-emitted x-rays. A comparison of the scale lengths of the Ha and x-ray 
emission along the minor axis confirms the more extended nature of the x-ray component: the Ha 
surface brightness along the minor axis is fit well by an exponential function, with a scale length 
of ~250 pc. For the most distant Ha emission (r ~ 1 kpc), this exponential can be approximated 
by a power law of slope —2, essentially the same power law exponent measured for the X-rays at 



a comparable radius ( Bregman, Schulman, fc Tomisaka 1995 ). Beyond this radius, the optical 



surface brightness falls more rapidly than does the x-ray surface brightness. 

Recent detailed modeling of the x-ray emission ( [Bregman, Schulman, fc Tomisaka 1995| ) 
suggests a temperature at large radii of only ~ 2 x 10^ K, implying an increasing role for thermal 
emission with radius. Similarly, hydrodynamic simulations find that the majority of the wind 



mass must be accumulated near the starburst region, not from evaporating halo clouds (Suchkov 



et al. 199^ ). It should be noted, however, that observations with the Ginga x-ray satellite have 



made the startling discovery of faint x-ray emission extending several tens of kiloparsecs from 
M82 ( [Tsuru et al. 1990| ). Hydrodynamic simulations have modeled this emission as shock-excited 
in nature, assuming that the outflow is much older, 5 x 10^ years (Tomisaka &: Bregman 1993). 



Although this observation has yet to be confirmed, the rapid radial decrease in the wind pressure 
and density could propagate the wind shock to large distances from its starburst origins. 

The true importance of shocks versus thermal emission can be estimated from optical line 
diagnostics, at least in the inner regions of the M82 outflow. As was pointed out in Section |3|, 
the flux- weighted [N II] /Ha ratio from the Fabry- Perot data is highly uniform and low in the 
inner kiloparsec of the outflow; values of 0.3-0.6 are typical. However, we must be careful to use 
line ratios for individual kinematic components when drawing conclusions regarding the physics 
of the gas, especially when the components have been modeled as distinct physical regions. The 
[N II] /Ha line ratio of the individual components reveals a similar low value across the spatial 
extent of split lines, except in the inner collimated zone, where a higher [N II] /Ha ratio is seen 
(~1.0), particularly in the low-velocity component. Although we were unable to resolve separate 
components in the [O III] observations, we note that the [O III] /Ha ratio exhibits a strong radial 
gradient, unlike the [N II] /Ha ratio. The [O III] /Ha ratio increases from a value of approximately 
0.03 at the center to 0.08 at a distance of ~750 pc. 

In order to investigate the importance of shock excitation for the optical filaments, we 
have compared the observed emission line ratios from the Fabry-Perot data with the standard 
emission-line galaxy diagnostic diagrams of [Veilleux &: Osterbrock 1987 . Although the small 



number of emission line diagnostics at our disposal limits our analysis, we can nevertheless make 
a rough assessment of the influence of shocks using the [O IIl]/H/3 versus [N II] /Ha diagnostic 
diagram (Fig. 14; from Fig. 1 of Veilleux fc Osterbrock 19871) . Using an H/3/Ha ratio of 0.25 for 



the outflow gas ( Heckman, Armus, &: Miley 199C| ), we see that the emission line ratios from the 



southern wind lobe of M82 rest in the region of the diagram for starburst galaxies, as expected. 
The ratios are comparable to those for most cooler H II regions and H II region models. This 
immediately suggests an emission mechanism such as photoionization for the filaments, particularly 
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near the nucleus, where the [O III] /Ha ratio is lower. As we move out in radius, however, shocks 
appear to become more important as an excitation mechanism, as the increasing [O III] /Ha ratio 



drives the locus in Figure 14 toward the region for non-thermally powered AGN. 



In order to more directly interpret our emission line fluxes in light of a shock mechanism, 
we have also compared the observed line ratios with a recent set of high-velocity shock models 
( Dopita &: Sutherland 19951 ). These models have been computed for shocks in the velocity range of 



150-500 km s ^; the deprojected gas velocity of the wind in M82 is estimated to be at the upper 
end of this range. Again using an H/3/Ha ratio of 0.25 for the outflow gas ( [Heckman, Armus, fc 



Miley 1990 ), a comparison with the [O IIl]/H/3 versus [N II] /Ha diagnostic diagram (Fig. 15; from 



Fig. 2b of Dopita &: Sutherland 1995|) shows that it is unlikely that the observed emission line flux 



from the inner outflow filaments arises entirely from shocks. There is simply not enough [O III] 
emission observed in the inner kiloparsec of the M82 outflow. However, the increasing [O III] /Ha 
ratio with distance from the nucleus suggests that shocks probably become important at larger 
radii, just as suggested by the observational diagnostic diagrams ([Veilleux fc Osterbrock T987D . 



Longslit optical observations have also reached the conclusion that the line ratios become more 
shock-like with increasing distance from the starburst (e.g., Heckman, Armus, fc Miley 1990| ), 



although this has often not included analysis of individual velocity components. 

Additional support for a photoionization mechanism for the inner optical filaments is provided 
by studies of the diffuse ionized medium (DIM) in NGC 891. In that galaxy, which has no outflow 
or other obvious sources of shock ionization, photoionization models have been used to understand 
the variation of line ratios with height above the disk plane. These models show that the [N II] /Ha 
ratio gradually decreases as the ionization parameter (the ratio of ionizing photons to gas density) 
increases. In contrast, the [O III] /Ha ratio should increase rapidly with ionization parameter 
( Sokolowski 199^ ). Regardless of the presence of shocks then, the low value of [N II] /Ha and the 



gradually increasing value of [O III] /Ha in the innermost filaments of M82 can be understood as 
the result of a gradual drop in filament density, relative to the number of ionizing photons from 
the starburst. In the outer filaments, however, the [N II] /Ha ratio begins to increase, presumably 
as a result of dilution of the radiation field in the expanding uncollimated bubble, as well as 
perhaps an increasing influence of shocks. 



Recent studies of the influence of halo dust on these line ratio trends in NGC 891 (Ferrara 



et al. 199^ ) point out that the [N II] /Ha ratio may appear artificially low near the disk due to 



dilution by scattered radiation from disk H II regions. This suggests that the low [N II] /Ha ratios 
for the inner filaments in M82 may be due in part to higher dust densities in the inner halo, 
scattering disk radiation from the nuclear starburst. This proposition corresponds well with the 
high levels of polarization detected from the filaments and the exponential nature of the observed 
halo, although the observed polarization levels in M82 (~ 10-15%) are much higher than those 
modeled in NGC 891 (~l-2%). 



Based upon these comparisons and our geometric models, we propose that the optical emission 
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from the inner kiloparsec of the M82 filament network is, at least partially, due to photoionization 
of the sides of the cavity created by the outflow. The hot gas in the wind itself would be quite 
transparent to the UV ionizing photons from the starburst region, allowing the entrained disk and 
halo gas to be illuminated directly. The tipped geometry of the outflow cones probably places the 
systemic side of each cone more directly in the path of the photoionizing radiation from the central 
starburst, explaining the higher fluxes seen in the low-velocity components of the wind. However, 
small regions of higher [N II] /Ha ratios in the individual velocity components suggest that a 
complex combination of shock and photoionization is probably required in the violent collimated 
zone, where the disk gas is being entrained and drawn upward by the hot wind just as it leaves 
the luminous starburst region. Although our small field of view and sensitivity limits restrict our 
analysis of the more extended optical filaments, we confirm a trend toward more shock-like line 
emission in the outer regions of the outflow. 



4-3.7. Comparison with UV Observations 

The inner regions of the outflow in M82 have been detected in UV images taken with the 
Ultraviolet Imaging Telescope (UIT; ptecher et al. 1992 ). Figure 16 compares our deep Ha image 
with the UV image observed in the 2490A band by the UIT. The southern outflow is clearly 
evident in ultraviolet emission, out to at least a kiloparsec. The northern outflow is only barely 
visible in the UIT data, presumably because of the high disk extinction at ultraviolet wavelengths. 
The projected morphology of the southern UV outflow is similar to that at Ha wavelengths, 
consisting of a relatively smooth, broad fan of emission. The ratio of the Ha and UV surface 
brightnesses (each in units of ergs cm~^ arcsec^^ A^^) in the outflow region is ~8.8 ( [Hennessy 



19961) . 



Previous analyses of the UIT data ( Maran et al. 1991 ; Hennessy 1993| ), as well as earlier 
balloon observations ( [Courvoisier et al. 1990| ; Blecha et al. 1990| ) have interpreted the minor-axis 
ultraviolet emission in M82 in terms of dust scattering of photons from the nuclear starburst by 
particles in the outflow. It has been suggested that filamentary structures can be seen in the 
UIT data aran et al. 1991| ), presumably because the higher relative densities in the filaments 
enhance scattering in those regions. The spatial distribution of the UV emission does not correlate 
well with the Ha, X-ray, radio continuum, or IR morphologies, although this may be due to 
variations in opacity ( Courvoisier et al. 19901 ). One difficulty with the scattering picture for the 
production of the southern UV outflow is the surprising non-detection of emission in the UIT 



1600A far- ultraviolet (FUV) band ( Hennessy 1993 ). The naive expectation for Rayleigh scattering 
would be to expect increased scattered flux at shorter wavelengths, even in the presence of 
substantial extinction along the line of sight. Ultraviolet observations of reflection nebula (e.g., 
Calzetti et al. 1995) and starburst galaxies (Calzetti, Kinney, & Storchi-Bergmann T99|) derive 



extinction curves that confirm this expectation. 



An alternative explanation for the UV light is "two photon" continuum emission from ionized 
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gas within the optical filaments. This emission arises from the spontaneous two-photon decay of 
the 2^5 level of H I and is commonly seen in AGN and regions of low-density ionized gas. The 
spectral distribution of two-photon emission is symmetric, in photons per frequency interval, about 
a peak at 2431 A (Osterbrock 1989), very close to the UIT 2490A band. Unlike the scattering 
function, two-photon flux decreases (~ A^^) toward shorter wavelengths (see also Hennessy 1996 ). 
Preliminary analysis of the [S II] AA6719,6731 lines in optical spectra of the outflow from the 
Keck telescope confirm previous studies (e.g., Rodriguez &: Chaisson 1980| ; Houck et al. 1984 ; 
Heckman, Armus, fc: Miley 1990 ) that the inner filament densities rarely exceed 500 cm~'^. As 



this is significantly below the critical density at which collisional effects influence the intensity of 



two-photon emission [ue ~ lO'^ cm-3; [Osterbrock 19891) , we may anticipate the filamentary nature 
of the outfiow in the UV imagery. 



The observed H/3/Hq; ratio of 0.25 for the outfiow gas (Heckman, Armus, & Miley 1990) 
implies an H/3/UV ratio of ~2.2. This is approximately a factor of ten greater than the ratio of 



H/3 to two-photon fiux predicted by photoionization models (~0.1; Dopita, Binette, & Schwartz 
1982; Dopita 19971 ), suggesting the presence of 2.5 mag more extinction at 2500A. The observed 



extinction law for starburst galaxies ( [Calzetti, Kinney, &: Storchi-Bergmann 1994 ; Fig. 5), 
combined with the observed Ha extinction factor of 1.8 in the outflow ( [Heckman, Armus, &: Miley 



1990) indicates that this level of obscuration is entirely reasonable. Shock- induced two-photon 
emission is almost certainly ruled out, as it would require extinction by a factor of ~ 4 to match the 
observed UV intensity, well above the value observed in the outflow regions. Spectral observations 
near the Balmer limit, ~3650A, where the two-photon continuum emission is enhanced relative to 
line emission, could be useful for determining the importance of this mechanism. 



5. Conclusions 

It is widely recognized that M82 is the prototype of galactic wind systems. Most models have 
concentrated on explaining measurements taken along isolated position angles close to the galaxy's 
minor axis. However, our new study has shown that the filamentary system associated with the 
outfiow is highly complex. Our simplest kinematic model requires at least two discrete structures 
for each side of the galaxy. The filaments are distributed in a network over these surfaces. More 
surprisingly, the axes of the outfiows on either side of the disk are aligned neither with each other 
nor with the spin axis of the disk. The gas excitation suggests that stellar ionization dominates 
close to the disk with increased dominance of shock excitation further along the outflow. The 
clearest signature of the biconic outflow geometry is seen in the line ratio maps which suggest that 
radiation is escaping from the disk along a channel excavated by the hot rarefied wind. We find 
evidence for a smooth line-emitting halo which we associate with the linearly polarized halo seen 
in broadband studies. We have observed a warm ionized medium throughout the inclined spiral 
disk. 

In our view, the most pressing issue is an explanation of the emission-line spectropolarimetry. 
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We strongly urge that these observations are repeated to deeper levels, at least to the point where 
the diffuse line-emitting halo is detected. Observations which only detect the bright filaments 
cannot remove the contribution from the halo along the line of sight. 



The disk-halo interaction in galaxies is a fundamental topic in its own rig ht (pioemen 1990|). 



We suggest that M82 may provide the best observational constraints on fountain flows (e.g., 
[Shapiro fc Benjamin 1991 ) once the origin of the filaments is fully understood. What happens 



to the entrained material that is lifted into the galactic halo? What is the relationship between 
the relativistic electron halo, the wind filaments, and the diffuse line-emitting halo? We suspect 
that answers can only come from future observations of x-ray, radio and millimeter emission at a 
resolution and sensitivity comparable to the present study. 
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Fig. 1. — Flux maps of M82 in the emission lines of (a.) Ha 6563A and {b.) [N II] 6583A, in units of 
ergs cm~^ sec~^, log-scaled between —15.5 and —13.0. The upper maps represent the total flux in 
the line at each pixel, while the lower panels represent the flux in the high- (HVC) and low- (LVC) 
velocity components of the line, in the regions where the line is split (designated by an outline in 
the upper panel). Splitting of the [N II] emission line is only observed south of the disk. 

Fig. 2. — Flux map of M82 in the emission line of [O III] 5007A, in units of ergs cm~^ sec~^, log- 
scaled between —16.0 and —14.2. The [O III] emission line profiles were fit with single components 
throughout the field. A contour map of the Ha flux has been superimposed. 

Fig. 3. — Velocity map of M82 in the emission line of Ha 6563A, in units of km s^^. The upper 
map represents the total velocity of the Ha-emitting gas at each pixel. In regions of multiple 
velocity components, the component velocities have been combined with flux weighting to yield an 
"average" velocity. An ellipse centered on the 2.2 fim nucleus ( [Waller, Gurwell, fc Tamura 1992 ) 



has been superimposed on the image. The lower panels represent the velocities in the high- (HVC) 
and low- (LVC) velocity components of the line, in the regions where the line is split (designated by 
an outline in the upper panel). Note that the velocity range is much greater for these component 
panels than for the total velocity map. 

Fig. 4. — Flux ratio map of M82 for the ratio of the flux in the emission line of [N II] 6583A to the 
flux in the emission line of Ha 6563A, log-scaled between —0.8 and -1-0.3. A contour map of the 
Ha flux has been superimposed. 

Fig. 5. — Flux ratio map of M82 for the ratio of the flux in the emission line of [O III] 5007A to 
the flux in the emission line of Ha 6563A, log-scaled between —2.2 and —0.5. 

Fig. 6. — Rotation curves from the major axis of M82, in order of increasing wavelength. Literature 



sources include: a.~b. McKeith et al. 199S, c. this paper, d.~e. Castles, McKeith, & Greve 1991, /. 



McKeith et al. 1993| , g. Beck et al. 1978, h. Loiseau et al. 1990. A systemic velocity of 208.7 km s ^ 



has been subtracted from the Fabry-Perot data. 

Fig. 7. — A profile of the peak Ha emission along a band parallel to and approximately 725 pc south 
of the major axis of M82, illustrating the bright filaments and underlying diffuse halo component. 
The dashed line represents a cut through an exponential halo model. Tick marks on the spatial 
axis are spaced at 1' intervals. 

Fig. 8. — A deep Ha image of M82. Filamentary structure is observed along the minor axis of the 
galaxy across the entire 6' field. A contour map of the ROSAT HRI image has been superimposed; 
the contours represent 0.5 to 4.0 events per pixel, in thirteen log-spaced increments. The x-ray 
map has been smoothed with a Gaussian of FWHM ~12". Tick marks are spaced at 1' intervals. 

Fig. 9. — Mosaic of 16 monochromatic frames from the Ha Fabry-Perot data cube. The field of 
view is 3^5. The frames are separated by ~45 km s~^, where the frame marked with an asterisk 
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corresponds roughly to the systemic velocity of the galaxy. The circular region in the southeast 
corner is a masked ghost from the star in the southwest corner. 

Fig. 10. — The radial velocity profile from fits to the Ha line components along the axis of the 
outflow (position angle of 150°), from both the Fabry- Perot data and the longslit observations 
of McKeith et al. 1995| . A systemic velocity of 200km s^^ has been subtracted from all values. 



The Fabry-Perot points have been interpolated over a band 7" wide through the 2.2 ^la nucleus. 
Dual Ha components are resolved in both the north and south outflow lobes. The shaded region 
represents a velocity cut along the axis of a two-cone Monte-Carlo simulation. 

Fig. 11. — Y-Z cuts along the axes of two Monte-Carlo models of the southern outflow in M82: a. 
a single truncated bubble and h. a pair of inclined cones. Dimensions are given in parsecs; arrows 
denote the direction toward the observer. 

Fig. 12. — Comparison of the observed Ha emission and a two-cone Monte-Carlo simulation of the 
southern outflow in M82. Panel a compares the model with the spatial distribution of Ha flux 
(contour map) and the regions of split Ha lines (grayscale). Panels b-d are two-dimensional spectra 
extracted along the axes drawn in panel a, and illustrate the multiple velocity components. The 
spectrum in panel h has been extracted along the axis of the cones, while panels c and d have been 
extracted perpendicular to the cone axis, at distances of 650 and 910 pc from the nucleus. Panel h 
is oriented with North at the top; panels c and d are oriented with West at the top. 

Fig. 13. — The Ha flux map from the optical Fabry-Perot data set overlaid with a contour map of 
the ROSAT HRI x-ray image. Tick marks are spaced at 1' intervals. The Ha image is displayed 
logarithmically between 10~^^'^ and 10"^^''' ergs cm~^ sec~^, while the x-ray contours represent 
to 5 events per pixel, in 7 equally-spaced increments. The x-ray map has been smoothed with a 
Gaussian of FWHM ~2". 

Fig. 14. — Reddening-corrected [O III]/H/? versus [N II] /Ha for a selection of H II regions (circles), 
H II region models (lines), starburst galaxies (triangles), and AGN (filled polygons). The position 
of the outflow gas in the southern lobe of the M82 wind is represented by the shaded region. The 
arrows are oriented in the direction of increasing distance from the nucleus and represent a radial 
extent of approximately 1 kpc. (Plot adapted from Veilleux &: Osterbrock 1987| .) 



Fig. 15. — Reddening-corrected [O III]/H/? versus [N II] /Ha grids for a selection of high- velocity 
shock models. Observations of Seyfert galaxies are shown as open circles; LINERs are shown 
as filled circles. The grid labeled 'Shock Only' includes only the emission from the shock, while 
the grid labeled 'Shock + Precursor' includes the contribution of preshock ionization. Emission 
characteristics of shocks with velocities of 150-500 km s~^ are computed. The position of the 
outflow gas in the southern lobe of the M82 wind is represented by the shaded region. The arrows 
ares oriented in the direction of increasing distance from the nucleus and represents a radial extent 
of approximately 1 kpc. (Plot adapted from Dopita &: Sutherland 1995[ ) 
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Fig. 16. — A deep Ha image overlaid with a contour map of the 2490A ultraviolet emission observed 
by the UIT. Tick marks are spaced at 1' intervals. 



Table 1. Useful observational values for the galaxy M82. All values are 
taken from de Vaucouleurs et al. 1991, except where otherwise noted. 



Designations 


M82, NGC 3034, UGC 5322, ARP 337, 3C 231 


Position (2000.0) 


a = 09^ 55?^ 54!0; d = 69° 40' 57" 


Galaxy type 


10, Irr II 


Dimensions 


llf2 X 4f3 (10.6 X 4.1 kpc) 


Position angle 


65° 


Disk inclination 


81?5 


Radial velocity 


203 ± 4 km s-i 


Distance 


3.63 ± 0.34 Mpc ^ 


Mass 


2.7 X IQi" Mq ^ 


Total magnitude 


B = 9.30 ±0.09 


Galactic extinction 


Ab = 0.13 


Total colors 


B - V = 0.89 ±0.01 




U - B = 0.31 ±0.03 


Radio luminosity 


Lradio = 10^^ ergs s'l ^ 


Infrared luminosity 


LiR = lO^'' ergs s~^ ^ 


Ha luminosity 


Lna 2 X lO-^i ergs s'^ 


X-ray luminosity 


Lx-ray = 2 X lO^" CrgS S"! ^ 



Lynds fc Sandage 1963 



''Freedman et al. 1994, For comparison purposes, we will use the standard value 
in the literature of 3.25 ± 0.20 Mpc ( Tammann fc Sandage 19'6^ ) in this paper. 

'^IWatson, Stanger, fc Griffiths 1984 



'' iHeckman, Armus, fc Miley 1990 
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Table 2. Theoretical and observational characteristics of the Fabry-Perot instruments. 







Hq:+[N II] Observations 


[0 III] Observations 


Constant 




theoretical 


experimental 


theoretical 


experimental 


Telescope /Instrument : 


CFHT + HIFI 




UH88 + HIFI 




Dtel 


telescope diameter 


3.6 m 




2.24 m 




fcam 


camera focal length 


90 mm 




50 mm 




Detector: 




IfA CCD + VISACOM 


CCD 






pixel size 


30 /Ltm 




24 iim 




P" 


plate scale 


0'.'86 pix-i 




0'.'846 pix-i 




rip 


CCD size 


256x256 pix 




432x432 pix 




t 


exposure time 


450 sec 




480 sec 




Etiilon: 




CF11T3 




SVl 




AA 


free spectral range (FSR) 


85 A 


87.9 A 


61 A 


55 A 


Nr 


reflective finesse 


60 


45-50 


60 




n 


order of interference 


75 


72 


91 


91 


Ay 


velocity resolution 


65 km s~^ 




60 km s-i 




Observations 










SX 


sampling increment 


0.96 A 


0.99 A 




0.67 A 


Sv 


sampling increment 


44 km s~^ 


45 km s~^ 




40.5 km s"^ 



Values include 2x2 binning. 
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Table 3. Parameters of the two-cone Monte-Carlo simulation which best reproduces the observed 

kinematics of the southern outflow in M82. 



Inner cone Outer cone 



inner height 


zi [pc] 





350 


outer height 


Z2 [pc] 


350 


800 


minimum diameter 


Di [pc] 


375 


390 


maximum diameter 


D2 [pc] 


400 


590 


opening angle 


a 


5° 


25° 


inclination angle 


b 


5° 


15° 


position angle 


pa 


150°-165° 


165° 


radial \"olc)(:it \" 


r{r) [km s-^; ([pc]) 


.72.-. 4 


- o.i:-i;- 
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